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We theoretically investigate spin transfer between a system of quasiequilibrated Bose-Einstein 
condensed magnons in an insulator in direct contact with a conductor. While charge transfer is 
prohibited across the interface, spin transport arises from the exchange coupling between insulator 
and conductor spins. In normal insulator phase, spin transport is governed solely by the presence of 
thermal and spin- diffusive gradients; the presence of Bose-Einstein condensation (BEC), meanwhile, 
gives rise to a temperature-independent condensate spin current. Depending on the thermodynamic 
bias of the system, spin may flow in either direction across the interface, engendering the possibility 
of a dynamical phase transition of magnons. We discuss experimental feasibility of observing a 
BEC steady state (fomented by a spin Seebeck effect), which is contrasted to the more familiar 
spin-transfer induced classical instabilities. 



PACS numbers: 72.25.Mk,72.20.Pa,75.30.Ds,03.75.Kk 

Bose-Einstein condensation (BEC) has been observed 
in a growing number of physical systems including 
trapped ultracold atoms and molecules [1 , semiconduc- 
tor exciton polaritons [2 , and microcavity photons [3 . In 
magnetic insulators, a quasiequilibrated BEC of magnons 
was created at room temperature by parametric pumping 
[T , which is especially intriguing as it represents the pos- 
sibility of phase transitions in spintronic devices. In the 
case of short-lived bosonic excitations such as polaritons, 
photons, and magnons, the system needs to be optically 
pumped to exhibit spontaneous condensation [5 . 

In magnetic systems, Gilbert damping of magnons is 
known to increase upon the introduction of an adjacent 
conductor [6 : If the magnet is made to precess, conduc- 
tion electrons may carry away spin upon colliding with 
the interface separating conductor and insulator, tilting 
the insulator's magnetization toward its axis of preces- 
sion. Known as spin pumping, this magnetic relaxation 
process is reciprocal to spin-transfer torque [II [8], by 
which the angular momentum and energy can be pumped 
back into the magnetic region [9 . We consider here the 
consequences of these reciprocal interactions on an in- 
sulator with inhomogeneous spatial fluctuations in the 
magnetization, in particular a system of Bose-condensed 
magnons similar to that mentioned above. In this Letter, 
we construct rate equations for spin transfer between a 
magnetic insulator and adjacent normal metal, and solve 
for the time-dependent spin accumulation in the metal 
and the phase behavior of the insulator. The main text 
is supplemented with a discussion of the thermodynam- 
ics of spin transfer in our system and proposal of possi- 
ble methods by which to detect the predicted dynamical 
phase transition. 

Let us consider the insulating ferromagnet subjected 
to a magnetic field B in the positive z direction and at- 
tached to a metallic conductor, as sketched in Fig. [l] 
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FIG. 1. The magnetic moments of insulator (left) atoms are 
coupled to the itinerant electrons of an adjacent conductor 
(right); an electron scatters inelastically off the interface, flip- 
ping its spin and creating or annihilating a magnon in the 
insulator. While coupling across the interface requires some 
degree of overlap between electrons in the conductor and lo- 
calized electron orbitals in the insulator, a net electron tun- 
neling between the two subsystems is prohibited, so that only 
spin density is transferred. The magnetic field in the insula- 
tor, and hence static magnetization, point in the positive z di- 
rection; for a negative gyromagnetic ratio the static spin den- 
sity is therefore oriented in the —z direction, so that magnons 
carry spin -\-h. 



Electrons in the ferromagnetic insulator are localized 
(typically in deep d oi f orbitals) near atomics sites, pre- 
cluding charge transport. The corresponding magnetic 
moments constitute individual degrees of freedom, which 
give rise to collective spin-wave excitations. Meanwhile, 
(s-character) electrons in the metal are considered com- 
pletely delocalized and noninteracting. We shall hence- 
forth denote the ferromagnetic subsystem as "left" or L, 
and the metallic conductor subsystem as "right" or R. As 
a starting point, we treat them as uncoupled so that the 
electronic state of the entire system is |m) = \mL)^\mR) . 
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\mL) is an eigenket of the linearized (i.e., noninteracting 
magnon) left Hamiltonian Hl] in other words, it is an 
element of the Fock space of Holstein-Primakoff (HP) 
magnons, each indexed by the mode number q. The 
magnon spectrum eq is gapped [min(eq) = egg > 0] by 
the presence of the magnetic field or anisotropy. Mean- 
while, \mR) is an element of electron Fock space and 
represents an antisymmetrized product of single-particle 
states corresponding to quasiparticle Hamiltonian 
each indexed by orbital quantum number k and spin a. 

Itinerant electrons in the conductor are coupled across 
the insulator-conductor interface to the magnetic mo- 
ments of the insulator by a generic exchange interaction. 
We suppose that this interaction Vint can be phenomeno- 
logically written in terms of creation (annihilation) oper- 
ators (cq) for free HP magnons and creation (annihi- 
lation) operators aj^^ (^ka) for conduction electrons: 



Vint = ^ Vqkk'Cqaj^,^ak; + H.c. , 



(1) 



qkk' 



where a =t or | denote electron spin in the +z or —z 
directions, respectively. Information about scattering off 
of the static component of the insulator magnetization is 
entirely contained in the conduction electron wavefunc- 
tion ip\^a- (x) , which we consider to have a finite albeit ex- 
ponentially vanishing extension into the insulator; more 
specifically, ?/^kcr (x) are eigenstates of the total mean-field 
Hamiltonian, including the interaction just on the inside 
of the insulator between the evanescent conduction elec- 
tron tails and the static z component of the insulator 
magnetization. We approximate the static component 
of the magnetization as spatially uniform in what fol- 
lows. The effect on conduction electron scattering due to 
the rotating magnetization component in the xy plane, 
i.e., Eq. ([T]), which we consider small in comparison to 
the static component, is responsible for spin pumping [6] 
and spin-transfer torque [7, 8 and treated perturbatively 
below. 

The first term on the right-hand side of Eq. ([T]) de- 
scribes a magnon (carrying spin up K) annihilating in 
the insulator to create a spin-down hole/spin- up elec- 
tron pair in the conductor, while its Hermit ian conjugate 
(H.c.) corresponds to a reverse electron spin- flip scat- 
tering off the insulator-conductor interface to create a 
magnon. The scattering amplitude Fqkk' is assumed to 
be a full matrix element describing this process. Notice 
that while energy is exchanged in this interaction, mo- 
mentum is not generally conserved. Moreover, this is not 
the only means by which conduction electrons can ex- 
change energy with the magnetic insulator: One could, 
for example, write down an inelastic scattering term of 
the form ~ Cq/Cqaj^/^akcr that conserves magnon num- 
ber (and therefore preserves the spin of the scattering 
conduction electron), which physically corresponds to a 
deviation of the spin-conserving part of the Hamiltonian 



from its mean-field form. Since such a process does not 
contribute to the flow of spins across the interface, how- 
ever, it becomes irrelevant when temperatures are main- 
tained by thermal reservoirs. It should also be noted that 
the presence of shape anisotropy generally gives rise to 
elliptical magnons. The elliptical magnon operators 6q 
and 6^ are linear combinations of circular magnon oper- 
ators Cq and c^, so that Cq and no longer diagonalize 
While our detailed analysis in the following assumes 
circular magnons, a finite magnon eccentricity is not ex- 
pected to significantly alter our findings qualitatively. 

The total Hamiltonian can be expanded as Htot = 
Hl + + Vint + Ht + ^env, whcrc Ht is a thermalizing 
Hamiltonian that contains magnon-magnon interactions 
and conduction electron-electron interactions, while Henw 
describes interactions between magnons and conduction 
electrons with their environments: magnon-phonon cou- 
pling, electron-phonon coupling, etc. Here we consider 
dephasing effects significant enough that coherence be- 
tween the left and right subsystems is destroyed and the 
density matrix for the entire system is always in the form 
Ptot = pL ^ pR' We further assert, subject to sufficiently 
fast thermalization in respective subsystems, that 



Tr[Pi?a^kaa'k'] = riF (/3R(ek - Pa)) ^^kk'^^a 



Tr[/5LC, 



riB {hi^ci- Pl))S, 



(2) 



where npix) = (e^ + 1)"-^ and nsix) = (e^ — 1)"-^ are 
the (quasiequilibrium) Fermi-Dirac and Bose-Einstein 
distributions, respectively, and ek (cq) is the electron 
(magnon) spectrum. Because each subsystem maintains 
internal equilibrium, magnons obey Bose-Einstein statis- 
tics while conduction electrons are described by a Fermi- 
Dirac distribution. Information about the allotment of 
spin and energy between them is now contained in the 
inverse temperatures and Pr^ the chemical potential 
Per for conduction electrons with spin a, and the effective 
magnon chemical potential pL (which does not have to 
vanish in a pumped system). Note that pl < Cgg, where 
egs is the ground-state magnon energy; the magnons be- 
come Bose-Einstein condensed when jj^l = ^gs- 

It is straightforward to calculate the spin current (per 
interfacial area A) j flowing into the insulator from the 
conductor in terms of temperatures and chemical poten- 
tials to lowest order in Vint using Fermi's golden rule: 
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where the ground-state, jgs, and excited, jex, magnon 
contributions are functions of the magnon chemical po- 
tential pl^ electron spin accumulation A/i = /i^ — /i^, 
and their temperatures Tl and Tr. In the thermody- 
namic limit, the spin-current density jgs, describing the 
rate of flow of ground-state magnons into and out of 
the insulator, is proportional to the number of ground- 
state magnons A/gs per insulator volume Vl, ^gs = 



3 



iVgs(ML,TL)/FL: 



(4) 



Here, gn is the Fermi- level density of states of conduction 
electrons and 
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where ep is the Fermi energy (assumed to be much larger 
than egs and temperature) and Vr volume of the conduc- 
tor. Note that the current density jgs is only present in 
the thermodynamic limit in BEC phase, fiL = ^gg. On 
the other hand, the spin-current density jex (carrying 
spin transfer via the excited magnon states) is present in 
both normal and BEC phases and, after some manipula- 
tions, can be written as 

Jex=2^ / de\V,^{e)\^{/^^-e)glgL{e) 

X [riB {hie - Ml)) - (/3i?(e - A/i))] , (6) 

in terms of the energy-dependent density of magnon 
states ^L(e). The (relatively weakly) energy-dependent 
quantity 
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contains information about inelastic transition rates in- 
volving excited magnons. 

The dynamics of spin flow across the interface are 
therefore determined by the sum of the condensate cur- 
rent density jgs, which is determined by spin accumula- 
tion in the conductor and the ground-state magnon en- 
ergy egs (and thus the applied magnetic field), and the 
thermal current density jex, which depends on both tem- 
perature and spin-potential biases. Note that sufficiently 
large spin splitting A/i in the conductor could, in princi- 
ple, drive spin density into the insulator until the required 
density of magnons is attained and the system undergoes 
Bose-Einstein condensation. In a recent experiment by 
Sandweg et al spin pumping into a metal by mag- 
netic insulator is driven by the presence of parametrically 
excited magnons; in addition, a spin current between the 
metal and insulator arises from a thermal gradient as dis- 
cussed above. The authors of Ref. [10 made use of the 
inverse spin Hall effect, wherein spin diffusion along a 
metal strip produces detectable Hall signal. Reciprocally, 
an electric current could be used to generate spin accu- 
mulation on the surface of a metal via the spin Hall effect; 
this surface spin accumulation may then drive magnons 
into the insulator [TT]. 
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FIG. 2. Behavior of predicted by the rate equation, 

^gs = jtot/hdL = jc/hdL —arigs/h. If jc had the sign opposite 
to that shown in the figure, the crossing point jc/o^dL would 
fall in the normal phase (ngs = 0), thus precluding a BEC 
formation. 



We henceforth focus on the regime where the temper- 
atures of both the left and right subsystems are fixed so 
that any energy gain or loss, independent of spin gain 
or loss, is completely absorbed or resupplied by thermal 
reservoirs. At fixed Tl the density of excited magnons 
Hex becomes a monotonic function of jul ^ ^gs alone. Let 
us further suppose that spin accumulation A/i in the right 
reservoir is independent of spin diffusion from the insu- 
lator and fixed. If the total density of magnons exceeds 
the critical BEC density Uc (corresponding to jj^l = Cgg), 
nex reaches and remains pinned at this value, nc, and 
only rigs is free to vary. In BEC phase, then, the time 
dependence of ngs is given by 



ngs(t) 



Jc 

adjL 



ngs(O) 



Jc 

adjL 



-xtjn 



(8) 



where the excited magnon fiux jc — jexifJ^L ^gs) is time 
independent, as long as /i^ is anchored by the condensate 
at egs, a = 27r (cgs - ^lAo'ii/dL, and dL = Vl/A 
is the magnetic layer thickness. The behavior of the 
Bose-Einstein condensed system thus falls into one of four 
regimes, as depicted in Fig. [2] In the first, A/i > egs (so 
that a < 0) and ngs(O) > jc/<^<^L, %s grows exponen- 
tially until saturating at a value ~ Ms/jhb (where Ms 
is the magnetization of the ferromagnet and /i^ is the 
Bohr magneton). In this case, magnon- magnon interac- 
tions become important ultimately and the system must 
be treated more carefully here. This is a realization of 
the "swaser" (i.e., a spin- wave analog of a laser) put for- 
ward in Ref. and observed in the context most similar 
to ours (in a magnetic insulator YIG) in Ref. [11 . In the 
second regime, A/i > egs but ngs(O) < jc/(^dL (requiring 
jc < 0), ngs decreases towards zero, and the system en- 
ters normal phase. The last two regimes (corresponding 
to jc > and jc < 0), which are of more interest to us. 
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/3r {A/i - egs) 

FIG. 3. When A/i < Cgs, the steady-state phase is insensitive 
to the initial condition for rigs, but depends on the tempera- 
ture bias Tl —Tr and the difference A/i — Cgs. As the sphtting 
A/i increases, the critical temperature for Tl increases until 
it equals Tr. Examples of time dependence in the normal 
and BEC phase regions are shown in the upper and lower left 
panels, respectively. When A/i > egs, depending on the ini- 
tial condition, the driven magnon system is either unstable or 
relaxes towards the normal phase. 

occur when spin splitting in the conductor is sufficiently 
small that A/i < egg and thus a > 0, as depicted in 
Fig. |3] Here, the steady-state phase no longer depends 
on the initial condition: When jc > 0, the magnons will 
Bose-Einstein condense (lower half of the main panel in 
Fig. |3|, and if jc < 0, normal phase with ngg = must 
eventually be reached (upper half of the main panel in 
Fig.[3]). 

In the normal phase (nex < ^c), Ml acquires time de- 
pendence, and the rate of change of the total number 
of magnons is ntot = ^ex = jex{t)/hdL- To illustrate 
these dynamics in a specific example, we consider a sim- 
ple model where the density of magnon states per unit 
insulator volume Vl has the form ^L(e) = SL(e/egs — 1)^ 
(with w > and Gl positive real number). In terms of 
the polylogarithm function 

Li^+i {z) = / dx^^:^ — - , (9) 

the density of excited magnons becomes 

ne^ = GLkBTL ^ , (10) 

where zl^Pl^^iil) = e^^'^^^"^^^^ is the effective magnon 
fugacity (with zl = I corresponding to a BEC). Assum- 
ing for simplicity that Vex{e) is energy independent and 
equal to Vg^^ one can see from Eq. (|6| that jex contains 
four terms, proportional to Ijiw+ii^L) / {PL^gs)^~^^ oc ngx, 
Li^+i(^i?)/(/3i^6gs)-+\ and Li^+2(^L/i?)/(/3L/i?egs)-+^ 



where zr = zl{Pl pR^fJ^L ^m)- At low tempera- 
tures {Pl^Pr ^ e~^) and assuming A/i < egg, the last 
two terms may be neglected in favor of the first two, 
leaving for an excited spin current: 

jex = Q^C^l(7 - ^ex) , (H) 

where 7 = nex(/^L ^ /^i?,ML ^ A/i) > 0. This yields 
time-dependent solutions to the excited magnon density: 

nex(t)=7+[^ex(0)-7]e-"'/^ (12) 

provided nex < ^c- Since A/i < egg, a > 0, and nex 
is decaying towards 7, irrespective of its initial condi- 
tion. If 7 < nc, thus, the insulator always remains in 
normal phase; when 7 > nc, on the other hand, the 
magnons eventually Bose-Einstein condense, and the sys- 
tem is henceforth described by Eq. (|8|. Notice that the 
conditions 7 ^ nc are (in the spirit of the aforementioned 
approximation) equivalent to jc ^ 0, which are com- 
pletely consistent with the conditions considered above 
for the system to settle in either normal or BEC phase 
as t ^ 00. 

If the insulator temperature Tl is left floating, the en- 
ergy flow between the two subsystems would give rise to 
the dynamics of Tl (supposing for simplicity Tr is still 
fixed). In the most extreme case, the insulator is allowed 
to exchange energy only with the conductor (and only 
by the electron-magnon scattering discussed above, ne- 
glecting phonon heat transfer), so changes in Tl are dic- 
tated by the rate at which energy is transferred across 
the barrier along with spin. The coupled rate equa- 
tions for energy and spin transfer can then be solved to 
give time-dependent solutions to the temperature Tl and 
the ground and excited magnon densities, nex and ngg. 
While this program is beyond our scope here, we may 
expect a significantly more complex phase diagram, with 
hysteretic features sensitive to the initial conditions and 
reentrant phase behavior. 

We conclude that BEC phase can be established under 
a steady-state transport condition when the ferromagnet 
is necessarily colder than the normal metal (thus facil- 
itated by a spin Seebeck effect [12 ) and the spin accu- 
mulation A/i is slightly below the spin-transfer torque 
instability (A/i ~ egg), in our model. (See Fig. js]) We 
have treated Gilbert damping as relatively unimportant 
compared to spin pumping, which is valid only for suffi- 
ciently thin magnetic films (~ 10 nm in case of permalloy 
with Gilbert damping ~ 10~^ [6] or YIG with damping 
~ 10~^ J3], but could be considerably thicker in high- 
quality YIG films that have Gilbert damping at least an- 
other order of magnitude lower [11 ). Considering that 
the unstable region (A/i > egg) has already been realized 
in practice ^ in a Pt/YIG bilayer spin biased by the in- 
verse spin Hall effect, and the spin-caloritronic properties 
[12] are presently under intense experimental scrutiny in 
such composites [lQl[T4], the experimental observation of 
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current-induced BEC phase in Pt/YIG hybrids appears 
very feasible. YIG film thickness larger than the char- 
acteristic de Broglie wavelength of magnons 1 nm at 
room temperature using standard YIG parameters p°5] ) 
would justify a three-dimensional treatment of BEC. A 
< 1 /im-thick YIG film with Gilbert damping < 10~^ like 
that employed in Ref. [iT thus appears adequate to our 
ends. Implicit in our discussion is the assumption that 
the magnon gas is dilute and can therefore be treated 
as noninteracting, aside from thermalization effects. In 
reality, these interactions must be accounted for. In this 
treatment, spectral properties would be somewhat modi- 
fied, but the essential behavior of the system is expected 
to remain unchanged. The emergent magnon superfluid 
properties [16] due to their interactions are left for a fu- 
ture work. 
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THERMODYNAMICS OF SPIN TRANSPORT 



Eq. (|6| suggests that spin flow to and from excited magnon states vanishes when there is no thermal or spin gradient, 
i.e., when Pl =^ Pr and Aja = jiL- However, when either of these conditions is not met, jex 7^ and spin (as well as 
energy) is transported across the insulator /conductor interface. In a steady state (i.e., zero spin current), in normal 
phase with thermal bias, Pr — Pl 0, a spin chemical potential difference 6fi = Afi — /j^l develops to oppose it: 



6fi ^ {(3r - Pl) 



where /3 = {Pl + Pr) /2, /i = (/iL + A/i) /2, = defiBi^)^ and the thermodynamic biases are assumed to be small 
(i.e. pR- pL <^P,S/j.<^fl). 

On the other hand, the condensed spin current jgs is independent of both Tl and T^, and, provided egg > A/i, 
always carries spin away from the conductor, irrespective of the temperature gradient between the two systems. The 
explanation for this behavior can be understood as follows. Consider a single tunneling event involving the creation 
(destruction) of a ground-state magnon (AA/gg = ±1) and the corresponding creation of a down- (up-) spin electron- 
hole excitation in the conductor {ANn = — AA^gg), which we call process A (B) in Fig. |4j The entropy change in the 
insulator associated with either process vanishes when the magnons form a BEC, so that the entropy change of the 
whole system is just dSu^ which can be found by enforcing energy conservation: 



A^t< 



Tr 



(egg - A/i) ANr . 



Thus, process B (A) is favored {ANr ^ 0) for tunneling events involving ground-state magnons when egg ^ A/i, in 
agreement with Eq. Q. Put differently, if A/i = the phase space of the conductor is unaffected with either the 
introduction of an up-spin excitation or the introduction of a down-spin excitation. However, process A requires the 
conductor to surrender an energy quantum egg to the insulator, whereas process B means a net gain in energy for the 
conductor; the overall entropy gain in the conductor (and therefore the entire system) is thus greater for process B 
than A. The zero-temperature version of this explanation is presented graphically in Fig. |4j 




FIG. 4. A down electron may relax the ferromagnetic insulator, carrying away the excess energy away in a scattering state 
above the Fermi surface cf (process B). An incident up electron on the Fermi surface, however, cannot transfer up spin to the 
insulator magnetization (process A), since such an energy-preserving process would raise the energy of the magnet, lowering 
that of the electron and therefore landing it below the Fermi surface, which is Pauli blockaded. Process B therefore dominates, 
and the insulator magnetization relaxes towards the easy axis. 
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DETECTION OF PHASE TRANSITION 

The BEC-normal phase transition presents some of the most interesting physics of the system, yet as can be seen 
from Fig. [3j it is difficult to discern from the total magnetization of the insulator alone: Whereas for fixed Tl the 
density of excited magnons plateaus as z^;, ^ 1, the rate of change of the total number of magnons til = fiex + ^gs 
remains always continuous function of time. The transition can, however, be observed by Brillioun light scattering, 
wherein the scattered light intensity scales quadratically with the lateral junction size if the ground-state condensate 
is indeed coherent [4 . 

Alternatively, electron spin resonance (or, for that matter, any spectroscopic probe of a coherent microwave radia- 
tion) can provide clear evidence of the presence of quasiequilibrated Bose-Einstein condensation of magnons. Consider 
a test-particle electron at a fixed distance r from the magnetic insulator. Provided that the electron experiences the 
insulator as a single quantum magnetic moment m, one may neglect details involving spatial fluctuations of the mag- 
netization and allow the two systems to interact via dipole-dipole coupling; the Hamiltonian describing the interaction 
is therefore of the form: 

^d-d = ^ rhiTijCFj , 

where Tij is a tensor that depends on r and cf is the electron spin operator. Supposing the electron, subjected to a 
strong applied magnetic field in the z direction, begins in the state the probability that quantum fluctuations in 
the magnetization m spin flip the electron is, to lowest order in T^j, 

Pt^;(t) = f dt' f dt" J2 TijT,j, {m it') m, it")) (tl aj \l) ai aj, |t) e'"^(*'-*") , 
Jo Jo ^j^,^, 

where ujz is the electronic Larmor frequency in the applied magnetic field. Choosing our coordinate system to coincide 
with the eigenbasis of Tij and for simplicity asserting cylindrical symmetry around the z axis (so that T^x = Tyy = T^), 
the transition probability becomes 

Pt^; (t) = /* dt' f dt"Tl {rh- {t') m+ {t")) e^-«(*'-*") = tTlS-+{w,) , 

JO ^0 

where S ^{uj) = J (ite*^^(m~(t)m+(0)) ex A^gg is the spectral density of magnetic oscillations in a steady state. The 

transition rate is thus proportional to Agg, which scales linearly with the lateral dimensions of our junction in BEC 
phase and is size independent in normal phase. This simple treatment is pertinent to the case when the magnons 
condense at q = 0. Otherwise (as is the case in YIG, for example) one needs to come up with means to couple 
coherently to magnetic fluctuations at a finite q (perhaps using some form of grating). 



